Gravitational Particle Production and the Moduli Problem 
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A theory of gravitational production of light scalar particles during and after inflation is investigated. 
We show that in the most interesting cases where long-wavelength fluctuations of light scalar fields 
can be generated during inflation, these fluctuations rather than quantum fluctuations produced after 
inflation give the dominant contribution to particle production. In such cases a simple analytical 
theory of particle production can be developed. Application of our results to the theory of quantum 
creation of moduli fields demonstrates that if the moduli mass is smaller than the Hubble constant 
then these fields are copiously produced during inflation. This gives rise to the cosmological moduli 
problem even if there is no homogeneous component of the classical moduli field in the universe. To 
avoid this version of the moduli problem it is necessary for the Hubble constant H during the last 
stages of inflation and/or the reheating temperature Tr after inflation to be extremely small. 
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I. INTRODUCTION 

Recently there has been a renewal of interest in gravita- 
tional production of particles in an expanding universe. 
This was a subject of intensive study many years ago, 
see e.g. ij. However, with the invention of inflationary 
theory the issue of the production of particles due to 
gravitational effects became less urgent. Indeed, gravita- 
tional effects are especially important near the cosmologi- 
cal singularity, at the Planck time. But the density of the 
particles produced at that epoch becomes exponentially 
small due to inflation. New particles are produced only 
after the end of inflation when the energy density is much 
smaller than the Planck density. Production of particles 
due to gravitational effects at that stage is typically very 
inefficient. 

There are a few exceptions to this rule that have mo- 
tivated the recent interest in gravitational particle pro- 
duction. First of all, there are some models where the 
main mechanism of reheating during inflation is due to 
gravitational production. Even though this mechanism 
is very inefficient, in the absence of other mechanisms 
of reheating it may do the job. For example, one may 
consider the class of theories where the inflaton potential 
V(<f>) gradually decreases at large <j> and does not have 
any minima. In such theories the inflaton field <fi does 
not oscillate after inflation, so the standard mechanism 
of reheating does not work ^,|| . To emphasize this un- 
usual feature of such theories we call them non-oscillatory 
models, or simply NO models ||. Usually gravitational 
particle production in such models lead to dangerous cos- 
mological consequences, such as large isocurvature fluc- 
tuations and overproduction of gravitinos jj| . In order to 
overcome these problems, it was necessary to modify the 
NO models J3| and to use the non-gravitational mecha- 
nism of instant preheating for the description of particle 



production Q). 

There are some other cases where even very small but 
unavoidable gravitational particle production may lead 
either to useful or to catastrophic consequences ||-|9| . For 
example, it has recently been found that the production 
of gravitinos by the oscillating inflaton field is not sup- 
pressed by the small gravitational coupling. As a result, 
gravitinos can be copiously produced in the early uni- 
verse even if the reheating temperature always remains 
smaller than 10 s GeV Another important example 

is related to moduli production. 15 years ago Coughlan et 
al realized that string theory and supergravity give rise 
to a cosmological moduli problem associated with the 
existence of a large homogeneous classical moduli field 
in the early universe Jl0[] . Soon afterwards Goncharov, 
Linde and Vysotsky showed that quantum fluctuations 
of moduli fields produced at the last stages of inflation 
lead to the moduli problem even if initially there were no 
classical moduli fields 0. Thus the cosmological mod- 
uli problem may appear either because of the existence 
of a large long-living homogeneous classical moduli field 
or because of quantum production of excitations (parti- 
cles) of the moduli fields. In || it was pointed out that 
the problem of moduli production is especially difficult 
in the context of NO models, where moduli are produced 
as abundantly as usual particles. 

Recently the problem of moduli production in the early 
universe was studied by numerical methods in ||, with 
conclusions similar to those of Ref . |Q . As we are going to 
demonstrate, the main source of gravitational production 
of light moduli in inflationary cosmology is very simple, 
and one can study the theory of moduli production not 
only numerically but also analytically by the methods 
developed in This will allow us to generalize and 

considerably strengthen the results of Refs. (?]]|]. 

In particular, we will see that in the leading approx- 
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imation the problem of overproduction of light moduli 
particles is equivalent to the problem of large homoge- 
neous classical moduli fields 0. We will show that the 
ratio of the number density of light moduli produced dur- 
ing inflation to the entropy of the universe after reheating 
satisfies the inequality 



x > TrH q 



3mAf 2 



(1) 



Here m is the moduli mass, M p <~ 2.4 x 10 18 GeV is the 
reduced Planck mass, and Ho is the Hubble constant at 
the moment corresponding to the beginning of the last 
60 e-foldings before the end of inflation. 

In the simplest versions of inflationary theory with po- 
tentials M 2 (j) 2 /2 or A0 4 /4 one has H ~ 10 14 GeV. In 
such models our result implies that in order to satisfy 
the cosmological constraint — < 10~ 12 one needs to 
have an abnormally small reheating temperature Tr ^ 1 
GeV. Alternatively one may consider inflationary models 
where the Hubble constant at the end of inflation is very 
small. But we will argue that even this may not help, 
so one may need either to invoke thermal inflation or to 
use some other mechanisms which can make the moduli 



problem less severe, see e.g. |ll|Jl2|| . 

In the next section we outline the classical and quan- 
tum versions of the moduli problem and explain how each 
of them can arise in inflationary theory. In section III we 
describe the results of our numerical simulations of gravi- 
tational production of light scalar fields during and after 
preheating. In particular we verify our prediction that 
the dominant contribution to particle production comes 
from long-wavelength modes which are indistinguishable 
from homogeneous classical moduli fields. Finally in sec- 
tion IV we analytically compute the production of these 
long wavelength modes and derive Eq.(jl|). This section 
also contains our concluding discussion. 



II. MODULI PROBLEM 

String moduli couple to standard model fields only 
through Planck scale suppressed interactions. Their ef- 
fective potential is exactly flat in perturbation theory 
in the supersymmetric limit, but it may become curved 
due to nonperturbative effects or because of supersym- 
metry breaking. If these fields originally are far from 
the minimum of their effective potential, the energy of 
their oscillations will decrease in an expanding universe 
in the same way as the energy density of nonrelativistic 
matter, p m ~ a~ 3 (t). Meanwhile the energy density of 
relativistic plasma decreases as a~ 4 . Therefore the rela- 
tive contribution of moduli to the energy density of the 
universe may quickly become very significant. They are 
expected to decay after the stage of nucleosynthesis, vi- 
olating the standard nucleosynthesis predictions unless 
the initial amplitude of the moduli oscillations xo is suf- 
ficiently small. The constraints on the energy density of 



the moduli field p x and the number of moduli particles 
n x depend on details of the theory. The most stringent 
constraint appears because of the photodissociation and 
photoproduction of light elements by the decay products 
of the moduli fields. For m ~ 10 2 — 10 3 GeV one has 



^ < 1(T 12 

s 



10 



-15 



(2) 



see |13|,[l4| and references therein. In this paper we will 
use a conservative version of this constraint, — < 10~ . 

If the field x 1S a classical homogeneous oscillating 
scalar field, then this constraint applies to it as well if 
one defines the corresponding number density of nonrel- 
ativistic particles x by the following obvious relation: 



Px 
m 



mx 



(3) 



Let us first consider moduli x with a constant mass 
m ~ 10 2 — 10 3 GeV and assume that reheating of the 
universe occurs after the beginning of oscillations of the 
moduli. This is indeed the case if one takes into account 
that in order to avoid the gravitino problem one should 
have Tr < 10 8 GeV. We will also assume for definitcness 
that the minimum of the effective potential for the field 
X is at x = 0; one can always achieve this by an obvious 
redefinition of the field %• 

Independent of the choice of inflationary theory, at the 
end of inflation the main fraction of the energy density of 
the universe is concentrated in the energy of an oscillating 
scalar field <fi. Typically this is the same field which drives 
inflation, but in some models such as hybrid inflation this 
may not be the case. We will consider here the simplest 
(and most general) model where the effective potential of 
the field <fi after inflation is quadratic, 



(4) 



After inflation the field oscillates. If one keeps the 
notation <f> for the amplitude of oscillations of this field, 
then one can say that the energy density of this field is 
given by p(4>) = \<fi . 

To simplify our notation, we will take the scale factor 
at the end of inflation to be ao = 1. The amplitude of 
the oscillating field in the theory with the potential (^) 
changes as 



(5) 



The field x does not oscillate and almost does not 
change its magnitude until the moment t\ when H 2 (t) = 
becomes smaller than m 2 /3. At that time one has 



Px m2 xl 
P4> ~ 6H 2 (t)M 2 



2M 2 



(6) 



This ratio, which can also be obtained by a numerical 
investigation of oscillations of the moduli fields, does not 
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change until the time tn when reheating occurs because 
p x and Ptj, decrease in the same way: they are propor- 
tional to a~ 3 . 

At the moment of reheating one has p^tn) = 
tt 2 N(T)T^/30, and the entropy of produced particles 
s = 27r 2 iV(r)T|/45, where N(T) is the number of light 
degrees of freedom. This yields 



Px Xo T R 



ms 3mM 2 



(7) 



Usually one expects Tr^> to ~ 10 2 GeV. Then in or- 
der to have < 10~ 12 one would need xo "C 10 _6 M p . 
However, it is hard to imagine why the value of the mod- 
uli field at the end of inflation should be so small. If one 
takes xo ~ M p , which looks natural, then one violates 
the bound ^ < 10~ 12 by more than 12 orders of mag- 
nitude. This is the essence of the cosmological moduli 
problem [f[o| . 

In general, the situation is more complex. During the 
expansion of the universe the effective potential of the 
moduli acquires some corrections. In particular, quite 
often the effective mass of the moduli (the second deriva- 
tive of the effective potential) can be represented as 



c 2 H z 



(8) 



where c is some constant and H is the Hubble parameter 
jfJI . Higher derivatives of the effective potential may ac- 
quire corrections as well. This leads to a different version 
of the moduli problem discussed in 0] , see also . The 
position of the minimum of the effective potential of the 
moduli field in the early universe may occur at a large 
distance from the position of the minimum at present. 
This may fix the initial position of the field x ancl lead 
to its subsequent oscillations. 

A sim ple toy model illustrating this possibility was 
given in 



V 



im 2 X 2 + ytf 2 (x 



Xo) 



(9) 



At large H the minimum appears at x — Xo] at small H 
the minimum is at x — 0- Thus one would expect that 
initially the field should stay at xoj an d later, when H 
decreases, it should oscillate about x = with an initial 
amplitude approximately equal to xo- The only natural 
value for xo m supergravity is xo ~ M p . This may lead 
to a strong violation of the bound (j|) . 

A more detailed investigation of this situation has 
shown Jl^] that one should distinguish between three dif- 
ferent possibilities: c ^> 1, e ~ 1 and c«l. 

If c > O(10), the field x is trapped in the (moving) 
minimum of the effective potential, its oscillations have 
very small amplitudes, and the moduli problem does not 
appear at all p2[ . This is the simplest resolution of the 
problem, but it is not simple to find realistic models 
where the condition c > O(10) is satisfied. 



The most natural case is c ~ 1. It requires a complete 
study of the behavior of the effective potential in an ex- 
panding universe. There may exist some cases where the 
minimum of the effective potential does not move in this 
regime, but in general the effects of quantum creation of 
moduli in this scenario [QJ|] are subdominant with re- 
spect to the classical moduli problem discussed above 
[Q,^6|, so we will not discuss this regime in our paper. 

Here we will study the case c -C 1. In this case the 
effective mass of the moduli at H 3> to is always much 
smaller than H, so the field does not move towards its 
minimum, regardless of its position. Thus if there is any 
classical field xo it simply stays at its initial position until 
H becomes smaller than m, just as in the case considered 
above, and the resulting ratio ^f- is given by Eq. ((?]). 

The moduli problem in this scenario has two aspects. 
First of all, in order to avoid the classical moduli prob- 
lem one needs to explain why xo ~ 10 _6 ^/^-Mp, which 

is necessary (but not sufficient) to have n x /s < 10 -12 . 
Then one should study quantum creation of moduli in an 
expanding universe and check whether their contribution 
to n x violates the bound n x /s < 10~ 12 . This last aspect 
of the moduli problem was studied in . 

In inflationary cosmology these two contributions (the 
contributions to n x from the classical field x arL d from its 
quantum fluctuations) are almost indistinguishable. In- 
deed, the dominant contribution to the number of moduli 
produced in an expanding universe is given by the fluc- 
tuations of the moduli field produced during inflation. 
These fluctuations have exponentially large wavelengths 
and for all practical purposes they have the same con- 
sequenc es as a homogeneous classical field of amplitude 

To be more accurate, these fluctuations behave in the 
same way as the homogeneous classical field <p only if 
their wavelength is greater than H^ 1 . During inflation 
this condition is satisfied for all inflationary fluctuations, 
but after inflation the size of the horizon grows and even- 
tually becomes larger than the wavelength of some of the 
modes. Then these modes begin to oscillate and their am- 
plitude begins to decrease even if at that stage m < H . 
To take this effect into account one may simply exclude 
from consideration those modes whose wavelengths be- 
come smaller than prior to the moment t ~ to" 1 
when H drops down to m. It can be shown that in the 
context of our problem this is a relatively minor correc- 
tion, so we can use the simple estimate xo = \J (x 2 )- 

In order to evaluate this quantity we will assume that 
c <C 1 and m <C H during and after inflation. This re- 
duces the problem to the investigation of the production 
of massless (or nearly massless) particles during and after 
inflation. In the next section we study this issue and show 
that in the most interesting cases where inflationary long- 
wavelength fluctuations of a scalar field can be generated 
during inflation, they give the dominant contribution to 
particle production. This allows us to reduce a compli- 
cated problem of gravitational particle production to a 
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simple problem which can be easily solved analytically. 



III. GENERATION OF LIGHT PARTICLES 
FROM AND AFTER INFLATION 

In this section we will present the results of a numerical 
study of the gravitational creation of light scalar particles 
in the context of inflation. Consider a scalar field x with 
the potential 



V( X ) = 7> K ~ t R ) X 2 



(10) 



where R is the Ricci scalar. In a Friedmann universe 
R = — -|(aa + a 2 )- The scalar field operator can be 
represented in the form 



1 



(2tt)- 3 / 2 



d A k 



a k Xk(t)e lkx + alxt(t)e 



where the eigenmode functions Xk satisfy 

2 



Xk + 3-x'fc + 



m 2 - £i? 



Xk = 0, 



(11) 



(12) 



By introducing conformal time and field variables defined 
as T) = J ^ , fk = axk eq. ( ^2| ) can be simplified to 

K + ulf k = Q (13) 

where primes denote differentiation with respect to r] and 



<4 



(14) 



The growth of the scale factor is determined by the evo- 
lution of the inflaton field (j> with potential V(<f>). In 
conformal time 



a 
a 



8na 



(<pl-a 2 V(<t>)) 



a- 2 dV(<P) 
2 — cp + a — — — 
a ' 



= 0. 



(15) 



(16) 



For initial conditions for the modes in the first 
approximation one can use positive frequency vacuum 
fluctuations fk = ~j^ e ~ lkt ^ see e -S- @- However, when 



describing fluctuations produced at the last stages of a 
long inflationary period, one should begin with fluctu- 
ations which have been generated during the previous 
stage of inflation. For example, for massless scalar fields 



minimally coupled to gravity instead of fk = — ^=e 



one should use Hankel functions (lTj : 



1 p —ikt 

2k 



fk(t) 



ia(t) H 
kV2k 



e Ht | exp 



i k 



-at 



(17) 



where H is the Hubble constant at the beginning of cal- 
culations. To make the calculations even more accu- 
rate, one should take into account that long-wavelength 
perturbations were produced at early stages of inflation 
when H was greater than at the beginning of the calcu- 
lations. If the stage of inflation is very long, then the 
final results do not change much if instead of the Han- 
kel functions ( [l7|) one uses fk — ~^j: e ~ lkt ■ However, if 
the inflationary stage is short, then using the functions 
fk = -j=e~ %kt considerably underestimates the resulting 

value of (x 2 )- 

At late times the solutions to Eq. ( |l3| ) can be repre- 
sented in terms of WKB solutions as 



fk(v) 



Pk(v) A-i I V uj k dr] 
— — P J ~v~ — 6 J 



(18) 



where cckiv) an( l Pk(v) pl & y the role of coefficients of a 
Bogolyubov transformation. This form is often used to 
discuss particle production because the number density 
of particles in a given mode is given by nk — \Pk(v)\ 2 an d 
their energy density is LOk n k- As we will see, though, the 
main contribution to the number density of x particles 
at late times comes from long-wavelength modes which 
are far outside the horizon during reheating. As long 
as they remain outside the horizon these modes do not 
manifest particle-like behaviour, i.e. the mode functions 
do not oscillate. In this situation the coefficients a and j3 
have no clear physical meaning. We therefore present our 
results in terms of the mode amplitudes |/fc(?7)| 2 , which 
as we will show contain all the information relevant to 
number density and energy density at late times. 

At late times when — ~ H < m the long wavelength 
modes of x will be nonrelativistic and their number den- 
sity will simply be given by Eq. (^). Moreover the very 
long wavelength modes which are still outside the hori- 
zon at late times (e.g. at nucleosynthesis) will act like a 
classical homogeneous field whose amplitude is given by 



(x 2 



2n 2 a 2 



dkk 2 \f k 



(19) 



It is these very long wavelength modes which will dom- 
inate and therefore the quantity of interest for us is the 
amplitude of these fluctuations. 



In our calculations we assumed that 



= c 2 H 2 with 



c <C 1; the results shown are for c = but we also did 
the calculations with c = .01 and found that the results 
were independent of c in this range. 

Figure [l] shows the results of solving Eq. (O) for a 
model with the inflaton potential V(4>) — . These 
data were taken after ten oscillations of the inflaton field. 
The vertical axis shows k 2 \fk\ 2 as a function of the mo- 
mentum k. The momentum is shown in units of the Hub- 
ble constant at the end of inflation. 
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FIG. 1. Fluctuations vs. mode frequency at a late time. 
The lower plots show runs which started close to the end of 



inflation, with initial conditions fk = 



ing times range from 



•/2k 



The start- 



1.5M P (lowest curve) up through 



4>a = 2M P . The highest curve shows the Hankel function so- 
lutions given by eqs. ( pp| ) and ( pl| ) As we see, calculations 
with fk = -j= e~ lht produce the correct spectrum at large k 
but underestimate the level of quantum fluctuations at small 
k. 



The different plots represent runs with different starts 
of inflation, i.e. with different initial values of <f>. They 
all coincide in the ultraviolet part of the spectrum, but 
the runs which started towards the end of inflation show 
a significant suppression in the infrared. This shows that 
fluctuations produced during inflation are the primary 
source of the infrared modes, which in turn dominate the 
number density. 

The curve on top shows the Hankel function solutions 
(|l7|), which give 



IM 



a 2 H 2 
2fc 3 



(20) 



for de Sitter space, i.e. for a constant H. In the figure, we 
have corrected this expression by using for each mode the 
value of the Hubble constant at the moment when that 
mode crossed the horizon. For the \<fr model shown here 
the appropriate Hubble parameter for each mode can be 
approximated as 



H k = 



2vrA 
~3~ 



1 / k 
—In 

7T \H„ 



(21) 



where <j) e and H e are the values of the inflaton and the 
Hubble parameter respectively at the end of inflation. 

Note that if the Hankel function solutions (|l7|) are used 
as initial conditions for a numerical run then they do not 
change as the modes cross the horizon, so the upper curve 
of the plot can also be obtained from such a run. Rela- 
tive to this upper curve it's easy to see how the numerical 
runs show suppression in the infrared due to starting in- 
flation at late times and choosing the initial conditions in 
the form fk = ^7= e~ l , and suppression in the ultravio- 
let due to the end of inflation. The latter suppression is 



physically realistic. The infrared suppression should oc- 
cur at a wavelength corresponding to the Hubble radius 
at the beginning of inflation. 

The different regions of the graph illustrate effects 
which occurred at different times. During inflation long 
wavelength modes crossed the horizon at early times. 
Thus the far left portion of the plot shows the modes 
which crossed earliest. They have the highest amplitudes 
both because they were frozen in earliest and because 
the Hubble constant was higher at earlier times when 
they were produced. The lower plots don't show these 
high amplitudes because inflation began too late for these 
modes to cross outside the horizon and be amplified. Far- 
ther to the right the curve shows modes which were only 
slightly if at all amplified during inflation. The far right 
modes were produced during the fast rolling and oscil- 
latory stages. These modes are not frozen and can be 
described meaningfully in terms of a and f3 coefficients. 
The regularized expression 



IM 2 



1 



1/3, 



k\' + 



Re 



-2i J u) k d-q 



(22) 



shows why the amplitudes of these modes oscillate as a 
function of k. 

In short inflationary fluctuations are primarily respon- 
sible for producing infrared modes and post-inflationary 
effects account for ultraviolet modes, but it is the infrared 
modes that were outside the horizon at the end of infla- 
tion which dominate the number density at late times. 
The earlier inflation began the farther this distribution 
will extend into the infrared, and the long wavelength 
end of this spectrum will always give the greatest contri- 
bution to the number density of \ particles. 

Our numerical calculations are similar to those of 
Kuzmin and Tkachev H. However, they took a rather 
small initial value of the classical scalar field </>, which 
resulted in less than 60 e-folds of inflation. As initial 
conditions for the fluctuations they used fk — -^^ e ~ lkt . 
They pointed out that the results of such calculations 
can give only a lower bound on the number of \ particles 
produced during inflation. Consequently, similar calcu- 
lations performed in |)| could give only a lower bound on 
the number of moduli fields produced in the early uni- 
verse. 

Our goal is to find not a lower bound but the com- 
plete number of particles produced at the last stages of 
inflation in realistic inflationary models, where the total 
duration of inflation typically is much greater than 60 e- 
foldings. One result revealed by our calculations is that 
the effects of an arbitrarily long stage of inflation can 
be mimicked by the correct choice of "initial" conditions 
chosen for the modes Xk after inflation. Instead of using 
the Minkowski space fluctuations fk = ~~rj: e ~ lkt used in 
H as well as in our numerical calculations, one should 
use the de Sitter space solutions (jl7|), with H corrected 
to the value it had for each mode at horizon crossing. 
Using these modes at the end of inflation is equivalent to 
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running a simulation with a long stage of inflation. 

Our numerical calculations confirmed the result that 
we are going to derive analytically in the next section 



(see also 



The number of x particles (n^ 



*<x 2 » 



produced during the stage of inflation beginning at (f> = 
</>o in the simplest model M 2 <fr 2 /2 is proportional to (pg, 
whereas in the model A</) 4 /4 it is proportional to 4>o- Thus 
the total number of particles produced during inflation 
is extremely sensitive to the choice of initial conditions. 
If one considers 4>o corresponding to the beginning of the 
last 60 e-folds of inflation, the total number of particles 
produced at that stage appears to be much greater than 
the lower bound obtained in ||. As we will see, this will 
allow us to put a much stronger constraint on the moduli 
theories than the constraint obtained in [Hi. 



IV. LIGHT MODULI FROM INFLATION 

The numerical results obtained in the previous sec- 
tion confirm our expectation that in the most interesting 
cases where long-wavelength inflationary fluctuations of 
light scalar fields can be generated during inflation, they 
give the dominant contribution to particle production. 
In particular, in the case of c -C 1, m <C H most mod- 
uli field fluctuations are generated during inflation rather 
than during the post-inflationary stage. These fluctua- 
tions grow at the stage of inflation in the same way as if 
the moduli field x were massless |17fl: 



d(x 2 ) 
dt 



4^' 



(23) 



If the Hubble constant does not change during inflation, 
one obtains the well-known relation 



<x 2 } = 



H 3 t 



(24) 



However, in realistic inflationary models the Hubble con- 
stant changes in time, and fluctuations of the light fields 
X with m -C H behave in a more complicated way. 

As an example, let us consider the case studied in the 
last section. Here inflation is driven by a field <f> with an 
effective potential V(4>) = j(j) 4 at cj) > 0. This potential 
could be oscillatory or flat for <f) < 0. We consider a light 
scalar field \ which is not coupled to the inflaton field 0, 
and which is minimally coupled to gravity. 

The field <fi during inflation obeys the following equa- 
tion: 



Here 



3H(/> = -Xcf 



2 V 3 M v 



(25) 



(26) 



where </>o is the initial value of the inflaton field 
this case Eq. (E3h reads: 



d(x 2 



dt 96V3V M 3 



exp 



-12 



Mnt 



The result of integration at large t converges to 

<x 2 > = -MM 2 . 

u ' 2 \24arM* J 



(27) 



In 



(28) 



(29) 



This result agrees with the results of our numerical in- 
vestigation described in the previous section. 

[From the point of view of the moduli problem, these 
fluctuations lead to the same consequences as a classical 
scalar field x which is homogeneous on the scale if -1 and 
which has a typical amplitude 



Xo 



A 01 



2 24ttM2' 



(30) 



A similar result can be obtained in the model V((f>) 



^~4> 2 - In this case one has |jj 



- MpMt. 



The time-dependent Hubble parameter is given by 



H = 



M 



which yields 



Xo 



M4>1 

8nV3M 2 ' 



(31) 



(32) 



(33) 



As we see, the value of xo depends on the initial value 
of the field </>. This result has the following interpreta- 
tion. One may consider an inflationary domain of initial 
size H~ 1 ((f>o). This domain after inflation becomes expo- 
nentially large. For example, its size in the model with 
V{<j)) = ^-<j) 2 becomes @ 



))exp 



AM 2 



(34) 



In order to achieve 60 e-folds of inflation in this model 
one needs to take cf>o ~ 15M P . This implies that a typi- 
cal value of the (nearly) homogeneous scalar field x m a 
universe which experienced 60 e-folds of inflation in this 
model is given by 



These two equations yield the solution ml 



Xo 



v¥ 2 



5M. 



(35) 
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In realistic versions of this model one has M ~ 5 x 
1(T 6 M p ~ 10 13 GeV @. Substitution of this result 
into Eq. (0) gives 



2 x 10 



-10 



(36) 



This implies that the condition n x /s < 10~ 12 requires 
that the reheating temperature in this model should be 
at least two orders of magnitude smaller than m. For 
example, for m ~ 10 2 GeV one should have Tr < 1 GeV, 
which looks rather problematic. 

This result confirms the basic conclusion of Ref. |t]] 
that the usual models of inflation do not solve the mod- 
uli problem. Our result is similar to the result obtained 
in H by numerical methods, but it is approximately two 
orders of magnitude stronger. The reason for this differ- 
ence is that the authors of Ref. || used a much smaller 
value of (f>o in their numerical calculations. Consequently, 
they took into account only the particles produced at the 
very end of inflation, whereas the leading effect occurs at 
earlier stages of inflation, i.e. at larger <f>. 

In general one can get a simple estimate of xo = \J (x 2 ) 
by assuming that the universe expanded with a constant 
Hubble parameter H during the last 60 e-folds of in- 
flation. To make this estimate more accurate one should 
take the value of the Hubble constant not at the end of in- 
flation but approximately 60 e-foldings before it, at the 
time when the fluctuations on the scale of the present 
horizon were produced. The reason is that the largest 
contribution to the fluctuations is given by the time when 
the Hubble constant took its greatest values. Also, at 
that stage the rate of change of H was relatively slow, so 
the approximation H — Hq ~ const is reasonable. Thus 
one can write 



2tt 



2tt 



This gives 



% > T R H o 



3m M 2 



(37) 



(38) 



In the simplest versions of chaotic inflation with poten- 
tials A/ 2 2 /2 or A</> 4 /4 one has H ~ 10 14 GeV, which 
leads to the requirement Tr < 1 GeV. But this equation 
shows that there is another way to relax the problem of 
the gravitational moduli production: one may consider 
models of inflation with a very small value of Hq . For 
example, one may have n x /s ~ 10~ 12 for Tr ~ Hq ~ 10 7 
GeV. 

However, this condition is not sufficient to resolve the 
moduli problem; the situation is more complicated. First 
of all, it is very difficult to find inflationary models where 
inflation occurs only during 60 e-foldings. Typically it 
lasts much longer, and the fluctuations of the light mod- 
uli fields will be much greater. This is especially obvi- 
ous in the theory of eternal inflation where the ampli- 



tude of fluctuations of the light moduli fields can be- 
come indefinitely large jl7|]. In particular, if the con- 
dition to 2 ~ m 2 + c 2 H 2 with c <C 1 remains valid for 
X ^ M p , then one may expect the generation of moduli 
fields x > M p . This should initiate inflation driven by the 
light moduli [9. Then the situation would become even 
more problematic: we would need to find out how one 
could produce baryon asymmetry of the universe after 
the light moduli decay and how one could obtain density 
perturbations Sp/ p ~ 10~ 4 in such a low-scale inflation- 
ary model. 

One may expect that the region of values of x where its 
effective potential has small curvature m 2 <C H 2 may be 
limited, and may even depend on H. Then the existence 
of a long stage of inflation would push the fluctuations of 
the field x U P to the region where its effective potential 
becomes curved, and instead of our results for xo onc 
should substitute the largest value of x f° r which m 2 < 
H 2 . In such a situation one would have a mixture of 
problems which occur at c <C 1 and at c ~ 1. 

Finally, we should emphasize that all our worries about 
quantum creation of moduli appear only after one makes 
the assumption that for whatever reason the initial value 
of the classical field x i n the universe vanishes, i.e. that 
the classical version of the moduli problem has been re- 
solved. We do not see any justification for this assump- 
tion in theories where the mass of the moduli field in the 
early universe is much smaller than H . Indeed, in such 
theories the classical field x initially can take any value, 
and this value is not going to change until the moment 
t ~ H^ 1 ~ m . The main purpose of this paper was to 
demonstrate that even if one finds a solution to the light 
moduli problem at the classical level, the same problem 
will appear again because of quantum effects. 

This does not mean that the moduli problem is unsolv- 
able. One of the most interesting solutions is provided 
by thermal inflation pl[ . The Hubble constant during 
inflation in this scenario is very small, and the effects 
of moduli production there are rather insignificant. An- 
other possibility is that moduli are very heavy in the early 
universe, m 2 = m 2 + c 2 H 2 , with c > O(10), in which 
case the moduli problem does not appear [|12| . The main 
question is whether we really need to make the theory so 
complicated in order to avoid the cosmological problems 
associated with moduli. Is it possible to find a simpler 
solution? One of the main results of our investigation 
is to confirm the conclusion of Ref. R] that the simplest 
versions of inflationary theory do not help us to solve the 
moduli problem but rather aggravate it. 

In conclusion we would like to note that the methods 
developed in this paper apply not only to the theory of 
moduli production but to other problems as well. For ex- 
ample, one may study the theory of gravitational produc- 
tion of superheavy scalar particles after inflation p,|] . If 
these particles are minimally coupled to gravity and have 
mass m <C H during inflation, then one can use our Eqs. 
(EkJ), ([33]) to calculate the number of produced particles. 
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These equations imply that the final result will strongly 
depend on i-e. on the duration of inflation. If inflation 
occurs for more than 60 Hubble times, the production of 
particles with m <C H is much more efficient than was 
previously anticipated. As we just mentioned, if c/)q is 
large enough then the production of fluctuations of the 
field x may even lead to a new stage of inflation driven 
by the field x Hi- O n the other hand, if m is greater 
than the value of the Hubble constant at the very end of 
inflation, then quantum fluctuations are produced only 
at the early stages of inflation (when H > m). These 
fluctuations oscillate and decrease exponentially during 
the last stages of inflation. In such cases the final number 
of produced particles will not depend on the duration of 
inflation and can be unambiguously calculated. We hope 
to return to this question in a separate publication. 
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